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Abstract 

This paper is concerned with application of the classical Youla-Kucera parameterization to finding a set of 
linear coherent quantum controllers that stabilize a linear quantum plant. The plant and controller are assumed to 
represent open quantum harmonic oscillators modelled by linear quantum stochastic differential equations. The 
interconnections between the plant and the controller are assumed to be established through quantum bosonic 
fields. In this framework, conditions for the stabilization of a given linear quantum plant via linear coherent 
quantum feedback are addressed using a stable factorization approach. The class of stabilizing quantum controllers 
is parameterized in the frequency domain. Also, this approach is used in order to formulate coherent quantum 
weighted and ML control problems for linear quantum systems in the frequency domain. Finally, a projected 
gradient descent scheme is proposed to solve the coherent quantum weighted Mi control problem. 

I. INTRODUCTION 

Coherent quantum feedbaek eontrol builds on the idea of eonstrueting a feedbaek loop from the 
intereonneetion of quantum systems through field eoupling [13]. This teehnique avoids loss of quantum 
information in eonversion to elassieal signals whieh oeeurs during measurement, as a direet eonsequenee 
of the projeetion postulate in quantum meehanies [16]. The eoherent quantum eontrol approaeh aims 
at developing systematie methods to design measurement-free intereonneetions of quantum systems 
modelled by quantum stoehastie differential equations (QSDEs); see for example [11], [17], [21]. 
Quantum-optieal eomponents, sueh as optieal eavities, beam splitters and phase shifters, make it possible 
to implement quantum feedbaek systems governed by linear QSDEs [20], [14], [21], provided the latter 
represent open quantum harmonie oseillators [5], [7]. This important elass of linear QSDEs models the 
Heisenberg evolution of pairs of eonjugate operators in a multi-mode quantum harmonie oseillator that 
is eoupled to external bosonie fields. The eondition of physieal realizability (PR) of a linear QSDE 
as an open quantum harmonie oseillator is organised as a set of eonstraints on the eoeffieients of the 
QSDE [11] or, alternatively, on the quantum system transfer matrix [24] in the frequeney domain. 
Coherent quantum feedbaek eontrol problems, sueh as stabilization and robust eontroller design, are of 
partieular interest in linear quantum eontrol theory [11], [21]. These problems are amenable to transfer 
matrix design methods [27], [28], [8], [21]. There are elassieal approaehes to eontrol synthesis problems 
for linear multivariable systems based on the transfer matrix of the system [31]. In partieular, one of 
the important strategies in eontroller design for sueh systems is the stable faetorization approaeh [25]. 
The eentral idea of the faetorization approaeh is to represent the transfer matrix of a system as a 
ratio of stable rational matriees. This idea gives rise to a methodology whieh leads to the solution of 
several important eontrol problems; see [25]. One of the fundamental results used in the faetorization 
approaeh to elassieal eontrol is the parameterization of all stabilizing eontrollers, whieh is known as 
the Youla-Kucera parameterization. The Youla-Kucera parameterization was developed originally in the 
frequency domain for finite-dimensional linear time-invariant systems using transfer function methods, 
see [29], [30], and generalized to infinite-dimensional systems afterwards [4], [22], [25]. The state space 
representation of all stabilizing controllers has also been addressed for finite-dimensional, linear time- 
invariant [18] and time-varying [3] systems. Eurthermore, the approach was shown to be applicable 
to a class of nonlinear systems [9], [19], [1]. In the present paper, we employ the stable factorization 
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approach in order to develop a quantum counterpart of the classical Youla-Kucera parameterization 
for describing a set of linear coherent quantum controllers that stabilize a linear quantum system. 
In particular, we address the problem of coherent quantum stabilizability of a given linear quantum 
plant. The class of stabilizing controllers is parameterized in the frequency domain. This approach 
allows weighted M 2 and coherent quantum control problems to be formulated for linear quantum 
systems in the frequency domain. In this way, the weighted M 2 and control problems are reduced 
to constrained optimization problems with respect to the Youla-Kucera parameter with convex cost 
functionals. Moreover, these problems are organised as a constrained version of the model matching 
problem [6]. Finally, a projected gradient descent scheme is proposed to solve the weighted M 2 coherent 
quantum control problem in the frequency domain. 


II. NOTATION 

Vectors are assumed to be organised as columns unless specified otherwise, and the transpose 
acts on matrices with operator-valued entries as if the latter were scalars. For a vector a of operators 
a\,...^ar and a vector h of operators hi,...,hj, the commutator matrix [a,b^] := ab^ — {ba^Y is an 
(r X 5)-matrix whose (j,A:)th entry is the commutator [aj,b]^ := ajbk — bkaj of the operators aj and bk- 
Furthermore, := ((■)*)^ denotes the transpose of the entry-wise operator adjoint (■)*. When it is 
applied to complex matrices, (■)’*' reduces to the complex conjugate transpose (•)* := ((■))^- Also, Ir 

\lr 0 ^ 


denotes the identity matrix of order r, and Jr := 


0 -Ir 


is a signature matrix. The Frobenius inner 


product of real or complex matrices is denoted by (M,V)p := Tr(M*V) and generates the Frobenius 

IRi Ri] 


norm || ■ ||f- Matrices of the form 


Ri Ri 


are denoted by A(7?i,7?2). The imaginary unit is denoted by 


\/—1, and the (j,A:)th block of a matrix F is referred to as F/^. The notation 


’ A 

B 

C 

D 

+ D with 


refers to a 


variable 5 G C. The conjugate system transfer matrix (F(— 5 ))* is written as F ( 5 ). The Hardy space of 
(rational) transfer functions of type p = 2 ,00 is denoted by M^, (respectively, JJ^M’p). The symbol 0 is 
used for the tensor product of spaces. 


III. LINEAR QUANTUM STOCHASTIC SYSTEMS 

The open quantum systems under consideration are governed by linear QSDEs which model the 
dynamics of open quantum harmonic oscillators. The input-output maps of such systems can be described 
in the frequency domain by using the transfer function approach [27], [28]. We will now outline this 
framework which is used as a basis of the frequency domain synthesis approach to quantum control 
presented in this paper. 


A. Open Quantum Harmonic Oscillators 

Corresponding to a model of n independent quantum harmonic oscillators is a vector a of annihilation 
operators ai,...,art on Hilbert spaces The adjoint aj of the operator aj is referred to as 

the creation operator. The doubled-up vector a of the annihilation and creation operators satisfies the 
canonical commutation relations (CCRs) [16] 



[a, a'*'] [a, 

K,at] K,aT] 


— Jn-) 


a 



( 1 ) 


We consider a linear quantum system whose dynamic variables are linear combinations of the 
annihilation and creation operators, acting on the tensor product space := 0 ... 

a := Fia- 1 -F 2 a* = [Fi F 2 ] a, (2) 
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where F\ and F 2 are appropriately dimensioned complex matrices. The relations dU) and Q imply that 

[d, a^] = F[a,=: 0, 

where F := A(Fi,F 2 ) G in accordance with the doubled-up notation [8], and the complex 

Hermitian matrix 0 of order In is the generalized CCR matrix [24]. Now, consider an n-mode open 
quantum harmonic oscillator interacting with an external bosonic field defined on a Fock space [20]. 
The oscillator is assumed to be coupled to m independent external input bosonic fields acting on the 
tensor product space '.= di® ■ ■ - where denotes the Fock space associated with the yth 

input channel. The field annihilation operators which act on form a vector 

Their adjoints that is, the field creation operators, comprise a vector ^*(0- The 

field annihilation and creation operators are adapted to the Fock filtration and satisfy the ltd table 


The 


Q dt in terms of the corresponding doubled-up vector £F\-a{t) := 

linear QSDEs, derived from the joint evolution of the n-mode open quantum harmonic oscillator and 
the external bosonic fields in the Heisenberg picture, can be represented in the following form [24], [8]: 


Mn(0 

<(0 


dd(t) = Ad(t)dt-f fid^:^n(t), (3) 

A£^oyiX,{t) = Ca{t)At-\-. (4) 

Here, the first QSDE governs the plant dynamics, while the second QSDE describes the dynamics of 


the output fields in terms of the corresponding doubled-up vector £^oxjx{t) := 
and creation operators acting on the system-field composite space Also, t 

g ^ C2«x2m^ g C2mx2«^ ^ ^ ^2mx2m 


of annihilation 


^out (j) 

le matrices A G 


A B 


-i&H - \@L*JmL -@L*J,nA{S, 0) 

C D 


L A(5,0) 


(5) 

where H = H* = A{Hi,H 2 ) G is a Hermitian matrix which parameterizes the system Hamiltonian 

operator ^a^Ha, the matrix L = A(Li,L 2 ) G specifies the system-field coupling operators, and 

S G is the unitary scattering matrix. 

B. Open Quantum Harmonic Oscillator in the Frequency Domain and Physical Realizability 

The input-output map of the open quantum harmonic oscillator, governed by the linear QSDEs ([3]) 
and (01), is completely specified by the transfer matrix [27], [28], [21], [8] which is defined in the 
standard way as 


r(i) := 


’ A 

B 

C 

D 


( 6 ) 


where the matrices A,B,C,D are given by the (5,L,//)-parameterization dS]). In view of the specific 
structure of this parameterization, not every linear QSDE, or the system transfer matrix with an 
arbitrary quadruple {A,B,C,D), represents the dynamics of an open quantum harmonic oscillator. This 
fact is addressed in the form of PR conditions for the quadruple {A,B,C,D) to represent such an 
oscillator; see [11], [17] for more details. 

Eor the purpose of the present paper, it is convenient to take advantage of the frequency domain 
version of the PR constraints on the system transfer matrices. Eor its formulation, we will need an 
auxiliary notion. 

Definition 1: The matrix A and the state-space realization dU) are said to be spectrally generic if the 
spectrum <7(A) has no intersection with its mirror reflection about the imaginary axis in the complex 

plane: ^(Ajn ( —(7(A)) =0, that is, A -f-v fiO for all eigenvalues A,v G <7(A). □ 


3 

















In particular, a spectrally generic matrix A can not have purely imaginary eigenvalues. Also, we will 
use a speeial elass of eomplex matriees of order 2m: 

:= {A(S,0) : S e is unitary}. (7) 

Lemma 1: [24, Theorem 4, p. 2040] Suppose F is a square transfer matrix of order 2m with a 
speetrally generic minimal state-space realization ®. Then F represents an open quantum harmonie 
oseillator in the frequency domain if and only if 

F {s)JmT{s) = Jm (8) 

for all 5 G C, and the feedthrough matrix D = F(oo) belongs to the set in ©. □ 

A transfer funetion F, satisfying the eondition ®, is said to be (7m,7OT)-unitary; see, for example, 
[12] and referenees therein. Sinee we eonsider this property for square transfer matriees (in whieh ease, 
([8]) implies that |detF(/(o)| = 1, and henee, F(z(o) is nonsingular for any (O G M), (7m,7„j)-unitarity is 
equivalent to its dual form [12], [21]: 

F(5)7„F>)=7„. (9) 

In view of ® and ®, the feedthrough matrix D in Lemma [H inherits (7„j,7OT)-unitarity DJ^D* = 
D*JmD = Jm from the transfer funetion F by eontinuity. Moreover, for an arbitrary matrix S G 
the matrix D = A(5,0) is (7m,-/m)-unitary if and only if S is unitary. In what follows, the argument s of 
transfer funetions will often be omitted for brevity. 


IV. FEEDBACK INTERCONNECTION 


We will now eonsider a linear quantum plant and a linear quantum eontroller with square transfer 
matriees P and K, respeetively, eaeh representing an open quantum harmonie oseillator in the frequeney 
domain. By analogy with similar struetures in elassical control settings, we partition the veetors 
and £/pout of the plant input and output field annihilation operator proeesses in aeeordanee with Eig. [B 




£/r 



( 10 ) 


Here .0)-, denote the veetors of annihilation operators of the input and output fields 

of the elosed-loop system, and the input and output of the eontroller, whieh eorrespond to the 
elassieal referenee, output, observation and eontrol signals, respeetively. In order to bring the feedbaek 
interconneetion to a standard format, the plant transfer matrix P is modified to as 


^11 ^12 
^21 ^22 


Cl 

C2 


-Pll Pl^ 
IAl_ _Pi3_ 

Pll P 23 
Pa\ P 43 


P\2 Pl4 

/’ 32 _ P34 

P 22 P 24 
P 42 P 44 


A 

[Bi B3] ' [B2 B4] 

1_1 1 




1 

1 

1 

1-1 1 

(N 

'C2 

C4 



£>21 0 ■ 

0 D21 

1 

1 

D22 0 

0 D22 


Bi B 2 


Dn 

D 21 


Du 

D 22 


( 11 ) 


The intereonneetion of the modified plant and the eontroller is shown in Eig. |2] The elosed-loop transfer 
matrix between the exogenous inputs and outputs of interest ean be ealeulated through the lower linear 
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fractional transformation (LFT) of the modified plant and the eontroller in the frequeney domain [31], 

[ 8 ]: 

G = + ^21 (12) 

Note that, similarly to the elassieal ease, the intereonneetion in Fig. |2l provides a general framework 




Pin 









Fig. 1. This diagram depicts the way in which the original plant P is modified to by partitioning the vectors of the plant input 
and output field operators in Gil. This modified structure allows for the connection to another linear quantum system which acts as the 
controller. 





S^y 


Fig. 2. This diagram depicts the fully quantum closed-loop system which is the interconnection of the modified quantum plant and 
the quantum controller K. The effect of the environment on the closed-loop system is represented by s^r- 


for the feedback interconnection of two quantum systems, one of which acts as the plant and the other 
as the controller. This framework includes the conventional coherent quantum feedback interconnection 
shown in Fig. [3] In the latter figure, the exogenous inputs and outputs of the closed-loop system are 


r\ 

n 



Zl 

Z2 


Fig. 3. This diagram depicts the way in which the quantum system r 3 , the concatenation of Ti and the feedthroughs, is formed by 
grouping the exogenous inputs to and outputs from the closed-loop system. 


grouped together. Here, r 3 (which is the concatenation of Fi and the feedthroughs) and r 2 represent 
the transfer matrices and K of the modified quantum plant and the quantum controller, respectively. 
Also, note that the PR conditions are usually formulated for the case when the number of exogenous 
inputs to the closed-loop system is not less than the number of outputs of the controller [11], [17]. 

V. PARAMETERIZATIONS OF STABILIZING CONTROLLERS 

For the purposes of Section |Vll we will now briefly review the Youla-Kucera parameterization 
of classical stabilizing controllers together with related notions. The latter include stabilizability, 
detectability, internal stability, coprime factorizations and matrix fractional descriptions (MFDs). Despite 
the quantum control context, these notions will be used according to their standard definitions in classical 
linear control theory [31], [25]. 
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A. Stabilizability of Feedback Connections 

Consider the two-two bloek of the modified plant transfer matrix ^ in (fTTI) given by 


-^22:= 


P22 

P24 


A 

B2 

P42 

P44 


C2 

F >22 


( 13 ) 


The following lemma provides a neeessary and suffieient eondition for the internal stability of the 
feedbaek system in Fig. [2l 

Lemma 2: [31, Lemma 12.2, p. 294] Suppose {A,B 2 jC 2 ) in ([T3]) is stabilizable and deteetable. Then 
the elosed-loop system in Fig. is internally stable if and only if so is the system in Fig. |4| □ 





Fig. 4. Equivalent stabilization diagram. 


B. Stable Factorization 

Let the transfer funetion 3^22 in (fT3l) have the following eoprime faetorizations over 

^22=NM-^ =M-^N, (14) 

where the pairs {N,M) and {N,M) of transfer funetions in speeify the right and left faetorizations, 
respeetively. Then there exist U,V,U,V G MML whieh satisfy the Bezout identities: 

VM-UN = f MV-NU = L (15) 


The following lemma provides suffieient eonditions for the existenee of stable eoprime faetors for the 
system ^ 22 - 

Lemma 3: [31, p. 318] Suppose {A,B 2 ,C 2 ) in ([T3]) is stabilizable and deteetable. Then the coprime 
factorizations of ^22 over described by ([T4l) . ([T5l) . can be chosen so as 


M U 
N V 

y -u 

-N M 


A + B2F 

B2 

-L 

F 

I 

0 

C2 + D22F 

D22 

I 


A+LC2 

— {B2+LD22) 

L ■ 

F 

I 

0 

C2 

—D22 

I 


(16) 

(17) 


where F e ^nd L e are such that both matrices A -|- B 2 F and A -|- LC 2 are Hurwitz. 

Furthermore, the systems in ([T^ . ([T7]) satisfy the general Bezout identity 


' y 

-u' 

M 

U 


7 

o' 

_-N 

M 

N 

y 


0 

I 


( 18 ) 
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C. The Youla-Kucera Parameterization 


The following lemma, based on [31, Theorem 12.17, p. 316], applies the results on the Youla-Kucera 
parameterization in the frequency domain to the closed-loop system being considered. 

Lemma 4: Suppose the block ^^22 of tho modified plant transfer matrix in (fTTI) has the coprime 
factorizations over describedJ)y (IT4l) . Also, let the auxiliary transfer matrices U,V,U,V G 

in ([T?]) be chosen so that = V^^U, which is equivalent to (fTSl) . Then the set of all stabilizing 

controllers is parameterized by 

K={U + MQ){V+NQ)-^ 

= iV + QN)-\U + QM)=hFTiOy,Q), (19) 

where the common parameter Q e of these factorizations satisfies 


det(y + Ae)(oo) ^0, 


( 20 ) 


and Oy 


UV-^ 

y-1 


-V-^N 


is an auxiliary system. 


□ 


In what follows, the class of stabilizing controllers will be parameterized using MFDs. However, 
they can also be parameterized in the LFT framework due to the relationship between MFD and LFT 
representations [31, Lemmas 10.1 and 10.2, pp. 241-242]. 


VI. QUANTUM VERSION OF THE YOUEA-KUCERA PARAMETERIZATION 

We will now employ the material of Sections |III]-0 in order to describe stabilizing coherent quantum 
controllers in the frequency domain. The following lemma represents (7^,7^)-unitarity in terms of the 
Youla-Kucera parameter Q from ([T9l) . 

Lemma 5: Suppose the controller transfer matrix K is factorized according to ([T9]) . Then K is 
unitary if and only if the parameter Q G satisfies 

^ + Q~A + A~Q + Q~nQ = 0 (21) 


for all 5 G C, where 

^:=U~J^U-V~J^,V, (22) 

A:=M~J^U-N~J^V, (23) 

n : = M~J^M - N~J^N. (24) 

Eurthermore, under the condition (l20l) . the feedthrough matrix K{oo) is well-defined and inherits 
unitarity from K. 

Proof: The (7^,7^)-unitarity condition for the controller 

K{s)J^K{s)=J^, (25) 

which must be satisfied for all 5 G C, is representable in terms of the right factorization from (IT9]) as 

{{U+MQ)iy+NQ)-^) J^{U+MQ){V+NQ)-^ =J^. (26) 

By the properties of system conjugation, (|26]) is equivalent to {U 4-MQ) J^{U + MQ) = (V -|- 
^Q) 4-NQ). After regrouping the terms, the latter equality takes the form 

u~j^,u - y~7^y + Q~{M~ j^u- v~7^y) 

+ {U~J^,M-V~J^N)Q 
+ Q~{M~- N~J^,N)Q = 0. 
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This leads to (|2T]) . with 4>, A, IT given by (I22l) - (l24l) . The faet that eondition (l20l) makes the feedthrough 
matrix K{oo) well-defined follows direetly from (fT9l) . The (7^,7^)-unitarity of ^(oo) is established by 
taking the limit in (1251) as 5 —)■ 0 °. ■ 

The proof of Lemma [5] shows that the eonstraint (l2TI) on the Youla-Kucera parameter Q inherits its 
quadratic nature from (l25l) . However, (|2T]) becomes affine (over the field of reals) with respect to Q 
in a p^icular case when IT = 0. In view of (fT4l) . the transfer function IT in (l24l) is representable as 
n = M {J^ — and hence, it vanishes if the block ^22 of the modified plant is 

unitary. Since (7^,7^)-unitarity (l25l) and its equivalent dual form KJ^K =7^ (cf. ([8]) and dH)) impose 
the same constraints on the square transfer matrix K, a dual condition to the one described in Lemma [5] 
holds for the left factorization of the controller in (fT9l) . This leads to a dual constraint on Q, which 
corresponds to (|2TI) . with A, IT being replaced with their counterparts expressed in terms of N, M, 
U, V. 

Theorem 1: Suppose the block ,^^22 of the modified plant transfer matrix ^ in (fTTI) has the coprime 
factorizations over MM ’00 described by (fT4l) . Also, let the transfer matrices Uy.UM G in (fTSl) 

satisfy the general Bezout identity (fT^ . Then the set of all stabilizing (7^,7^)-unitary controllers K 
with a well-defined feedthrough matrix K{oo) is parameterized by (fT9l) . where the parameter Q belongs 
to the set 

^ := [Qe MMoo satisfying (f20j) and fl2Tjl }. (27) 

Proof: This theorem is proved by combining Lemmas |4] and [51 Indeed, since the underlying coprime 
factorizations are assumed to satisfy the general Bezout identity (flSl) . then (1211) can be applied to the 
common parameter Q in (fT9l) in order to describe all stabilizing (7^,7^)-unitary controllers K. Their 
feedthrough matrices K{oa) are well-defined provided the additional condition (l20l) is also satisfied. The 
resulting class of admissible Q is given by (|27]) . ■ 

Theorem [H provides a frequency domain parameterization of all stabilizing (7^,7^)-unitary controllers 
with a well-defined feedthrough matrix and leads to the following theorem. 

Theorem 2: Under the assumptions of Theorem [H the MFDs (fT9l) describe a set of stabilizing PR 
quantum controllers K, where the parameter Q belongs to the following class defined in terms of 
([27]) and ©: 

^ [Q ^ ^ : K in ([T9|) is spectrally generic, 

and 7f(oo) e }. (28) 

Proof: The assertion of the theorem is established by combining Theorem [T] with the frequency 
domain criterion of PR provided by Lemma [TJ ■ 

Theorem |2] parameterizes a subset of stabilizing PR quantum controllers in the frequency domain 
through the representation (fT9l) and the set £2 in (l28l) . This subset of controllers does not exhaust 
all stabilizing coherent quantum controllers. However, the discrepancy between these two classes of 
controllers is only caused by the technical condition of spectral genericity which comes from Lemma [B 

VH. COHERENT QUANTUM WEIGHTED Mi AND CONTROE PROBEEMS IN THE 

EREQUENCY DOMAIN 

The following lemma, which is given here for completeness, employs the factorization approach in 
order to obtain a more convenient representation of the closed-loop transfer function. 

Lemma 6: Under the assumptions of Theorem [T1 for any stabilizing controller K parameterized by 
([T9l) . the corresponding closed-loop transfer matrix G in (fT^ is representable as 

G = To + TiQT2, (29) 

where 

To:=Mn + Mi2UMM2u Ti:=MnM, T 2 —MM 21 . (30) 
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Lemma allows the following coherent quantum weighted and control problems to be 
formulated in the frequency domain. 

A. Coherent Quantum Weighted M 2 Control Problem 

Using the representation (l29l) . we formulate a coherent quantum weighted M 2 control problem as the 
constrained minimization problem 

£:=||l^outGlUi„||i=||To + TieT2||i^min (31) 

with respect to 2 G where the set ^ is given by (1281) . Here, 

To:=lUoutWm, Ti:=lLoutri, T 2 := TzlUin, (32) 

where Tq, Ti, T 2 are defined by (l30l) . Also, VLin, Wout G MM 00 are given strictly proper weighting transfer 
functions for the closed-loop system G which ensure that Tq-{-T\QT 2 G M 2 - The ,^-norm 11-112 
is associated with the inner product (ri,r 2 ) := ^(ri(z(o),r 2 (ftt)))pd( 0 . By using the standard 
properties of inner products in complex Hilbert spaces [23], the cost functional E in (|3TI) can be 
represented as 

£ = ||To||i + 2Re(To,2)-h(2,Ti2f2), (33) 

where 

fo:=TiToT2, fi:=TiTi, f2:=T2T2. (34) 

In comparison to the original coherent quantum LQG control problem [17], the coherent quantum 
weighted M 2 control problem (|3T]) allows for considering the cost of the unavoidable quantum noise 
fed through the plant by the controller. 

B. Coherent Quantum Weighted ML Control Problem 

Similarly to (l3T]) . a coherent quantum weighted ML control problem is formulated as the constrained 
minimization problem 

||G||^= ||To + Ti2T2||oc^min (35) 

with respect to 2 G where the set ^ is defined by (l28]) . Here, To, Ti and T 2 are given by (l3^ . 
where, this time, the weighting transfer functions VLin,VLout £ MM ’00 are not necessarily strictly proper. 
Recall that the norm in the Hardy space M’oo is defined by ||r||cx, supg,gRC7max(r(ztt))), where O'niax(-) 
denotes the largest singular value of a matrix. Note that both problems (ISTI) and (l35]) are organised as 
constrained versions of the model matching problem [6]. Since the M ’2 control problem is based on 
a Hilbert space norm, its solution can be approached by using a variational method in the frequency 
domain, which employs differentiation of the cost E with respect to the Youla-Kucera parameter Q and 
is qualitatively different from the state-space techniques of [26]. 

VIH. PROJECTED GRADIENT DESCENT SCHEME EOR THE COHERENT QUANTUM 

WEIGHTED M ’2 CONTROE PROBEEM 

Suppose the set L2 in (|2^ is nonempty, and hence, there exist stabilizing PR quantum controllers for 
a given quantum plant. By using the representation (|3^ and regarding the transfer function Q G MM ’00 
as an independent optimization variable, it follows that the first variation of the cost functional E in 
(|3T]) with respect to Q can be computed as 

5T = Re(VT,52), VT := 2(To-hfi2T2), (36) 

where use is also made of (l34l) . In order to yield a PR quantum controller, Q must satisfy the constraint 
(1271) whose variation leads to 

52~(A+n2) + (A~+2~n)52 = 0. (37) 
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In view of the uniqueness theorem for analytie functions [15, pp. 369-371], the resulting constrained 
optimization problem can be reduced to that for purely imaginary s = ico, with (O G M. The transfer 
matrices dQ, satisfying (1371) at frequencies CO from a given set C M, form a real subspace of transfer 
functions 

y:={Xe^je^: (X*(A+n<2) + (A*+<2*n)X)|.^ = 0}. (38) 

For practical purposes, the set Q. is used to “discretize” the common frequency range of the given 
weighting transfer functions VFin, Wout in the coherent quantum weighted control problem (l3TI) . A 
numerical solution of this problem can be implemented in the form of the following projected gradient 
descent scheme for finding a critical point of the cost functional E with respect to Q subject to (l2TI) at 
a finite set of frequencies Q: 

1 ) initialize Q G so as to satisfy (l2TI) . which yields a stabilizing PR quantum controller; 

2) calculate 'VE{ico) according to (l3^ for each frequency (0 G Q; 

3) compute 5Q{i(o) = —a?vo]y{VE{i(o)) by using a projection onto the set and a parameter 

a > 0; 

4) update Q io Q-\- 5Q, and go to the second step. 

The gradient projection Proj^(V£') onto the set in (IMI) is computed in the third step of the algorithm 
by solving a convex optimization problem on a Hilbert space with the direct sum of the Frobenius inner 
products of the projection errors at frequencies ft) G This computation, which will be discussed 
elsewhere, also involves the interpolation of transfer functions; see [2], [10] for more details. The 
discrete frequency set Q. and the step-size parameter a can be chosen adaptively at each iteration of 
the algorithm. The outcome of the algorithm is considered to be acceptable if Q belongs to the set ^ 
defined by (|2^ of Theorem [2| In the case when Q satisfies all the conditions in (|28]) except for the 
spectral genericity, slightly different weighting matrices can be used in order to remedy the situation. 

IX. CONCLUSION 

The set of stabilizing linear coherent quantum controllers for a given linear quantum plant has been 
parameterized using a Youla-Kucera factorization approach. This approach has provided a formulation of 
coherent quantum weighted M 2 and ML control problems for linear quantum systems in the frequency 
domain. These problems resemble constrained versions of the classical model matching problem. A 
projected gradient descent scheme has been outlined for numerical solution of the coherent quantum 
weighted M 2 control problem in the frequency domain. The proposed framework can also be used to 
develop tractable conditions for the existence of stabilizing quantum controllers for a given quantum 
plant, which remains an open problem. This is a subject of future research and will be considered in 
subsequent publications. 
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